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DYNAMICS AND IDENTIFICATION OF FLEXIBLE AIRCRAFT 


William R. Wells 
University of Cincinnati 


INTRODUCTION 


The practical necessity of removing all excessive structural 
weight through either conventional design practices or an active 
aerodynamic control system, has resulted in vehicles that are more 
aeroelastic. However, most existing parameter estimation methods 
used in the study of the stability and control and handling proper- 
ties of the aircraft stop short of explicity identifying important 
aeroelastic parameters that affect the aircraft dynamics. Previously, 
structural motion was considered as part of the measurement noise and 
filtered rather than modeled dynamically. 

The purpose of this paper is twofold: (1) to present, in detail, 

a development of the equations of motion for an elastic aircraft in- 
dicating, when appropriate, the difference and similarities with the 
dynamics of the rigid airplane, (2) discuss the additional computational 
difficulties due to the inclusion of added stability derivatives 
arising from consideration of aeroelastic effects. 

EQUATIONS OF MOTION 


Rigid Body Motion 


The first set of equations to be obtained are the "rigid body 



equations of motion" for the elastic airplane. These result from the 
application of Newton's laws of motion to the aircraft as a whole and 
describe the motion of the aircraft center of mass relative to an 
inertial system. 

Consider the elemental mass and surface area shown in Figure 1. 
The externally applied forces acting on these are a body force, R, 
and a surface force F. The body force will be assiimed to be gravi- 
tational and the surface force aerodynamic. The thrust and control 
surface forces can be added independently after the equations are 
formulated under the assumption that these quantities do not contri- 
bute to the elastic deformation of the structure. In this figure, r 
is the location of the mass element relative to the center of mass 
before the aircraft distorts and d represents the elastic displacement 
of the particle due to distortion. 

Conservation of linear and angular momentum dictate that the 
following equations be satisfied: 

J p dV = J R dV + / F dS (1) 

V VS 

/ pr ' X dV = / r ' X R dv + J r ' X F dS (2) 

V VS 

where V cind S refer to volume and surface quantities respectively. 

The time rate of change of the vectors in the left-hand sides of 
equations (1) - (2) are relative to the inertial earth fixed system. 


2 





system 

system 


Figure 1. - Definition of axis systems used to describe the aircraft 
deformation. 
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If the vector r' is replaced by r^ + r in equation (1) , the 
result is 


d 

dt 


dr' 

/•p dV + 


d [ - 

dt J P 


dr 

dV = / R dV + / F dS 
V S 


The center of mass is located by the vector 



dV 


From the definition of the center of mass location, it follows that 


/ p r dV = 0 
V 


and 


d 

dt 


/ p r dV 
V 


~ dr 

/ P at av 

V 


These results allow equation (1) , 
expressed as 
dV 

M = / FdS+Mg. 

s 


the linear momentum equation, to be 


(3) 


where 


V H 


dt 


R = p £ 
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similarly, equation (2) can be written in the form 


£6 



M2. - / F dS] 
S 


+ 


d 

dt 


/ p r X 
V 


dr 

dV = / r X F dS 
dt s ~ 


The bracketed term in this equation vanishes due to equation (3) ; 
consequently, the conservation of angular momentvun equation reduces 
to 




dr 

3 TT dV = / r X F dS 
dt s “ 


(4) 


At this point in the development it is of interest to note how 
coupling between the elastic deformation and rigid body motions occurs. 
In equations (3) -(4) both vectors F and r depend on the shape of the 
aircraft and, hence, upon the elastic deformation represented by d. An 
objective of the analysis is to minimize, as much as possible, this 
coupling without sacrificing accuracy. 

Next, equations (3) -(4) will be written in terms of velocities 
and accelerations relative to the body axis system which translates 
and rotates relative to the inertial system with angular velocity 
The reason for doing this is to illustrate the similarity of this 
overall motion or axis system motion to the motion of the rigid air- 
craft and to indicate how elastic effects are treated. 

If and { ) represent the rate operators as measured in in- 
ertial and body axis systems, respectively, then 
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= V + to X V 
at — c — — c 

dr 

= r + u X r 

Further, let 

u = Pi + + Rk 

^ = Ui + V2 + Wk 

2 = -gsine^ + 9 coseaincj)^ + g cos6cos(^k 

/ FdS=Fi+Fi+Fk 
^ — X- y-^ z- 

/ r X F dS = M„i + M, + M„k 
Q — — X— z— 

where i, j_, and k are unit base vectors in the body fixed axis system. 

In scalar form, the linear momentum equations (Equation (3)) are 
written as 

MU + M(QW-RV) = -Mgsine + F^ 

MV + M(RU-PW) = Mgcosesinf +F^ 

MW + M(PV-QU) = Mgcosecosifi + F^ (5) 

It can be seen that in these equations, the rigid body motion and elastic de- 
formations are dynamically uncoupled; that is, the effects of elastic defor- 
mation enter the equations only in the right-hand side through the 
applied force terms. 

The angular momentum equation, from equation 4, can be rewritten 
as 
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(6) 


^^/Mrxd + rx (uxr)] dV = M 
V 

• • - 

Where the identities r=d and M = j r x F dS have been used. Note 

S 

that the rigid body motion and elastic deformation are dynamically coupled in 
this equation due to the presence of r and d in the left hand side 
of the equation. 

The dynamic coupling expressed in equation (6) can be eliminated 

by the adoption of two approximations; the rate of elastic displacement 

is slow relative to the velocity ^ x r and the total angular momentum 

about the center of mass of the deformed aircraft can be computed with 

r = r as if the aircraft is undeformed. 

The first approximation assumes that /prxddV-0 which is 

V 

also an approximation consistent with the normal mode method to be 
discussed later. This allows the angular momentum equation to be 
written as 


/ p r X (o) X r) dV = M (7) 

dt ^ - - - 

The second approximation, £ = £. allows the angular momentum 
equations to be written in terms of moments of inertia of the j ig 
shape and allows equation (7) to be written as 

^ / p r x (u X r) dV = M (8) 

at ^ 

For a constant mass, symmetric aircraft with the x y z axis as 
principal axes. Equation (8) can be expressed in scalar form as 
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p 


I 


XX 


I 


yy 


Q 


1 

zz 


R 


-I R-I PQ+(I -I )RQ = M 
xz xz ‘ zz yy ^ X 

+ (I -I )PR+I (P^ - R^) = M 
XX zz xz ' ' y 

-I P+(I -I )PQ + I QR=M 
XZ yy XX ^ xz “ z 


(9) 


In either equation (8) or (9) it is noted that the effects of 
elasticity enter through the applied moments. 

Equations (5) and (9) can be expressed in terms of perturbed 
state variables u, r, w, p, q, etc. from a reference flight condition 
^1' '^l' '^l' ^1' ^1' sake of brevity, the perturbed form of 

equations (10) and (11) are presented for the special reference con- 
ditions of zero roll (P^^ = 0) and level flight (w^^ = 0) as 


Mu + M(Q^w - R.v - V r) + Mge cose, = ^ 

± -L ± P 1 X 

Mv + M(Rj^u + Vj^r) + Mg (^sin0j^sin(t)j^ -^p cos9^cosi(ij^) = f^ (10) 
Mw + M(V,p - Q,u - U,q) + Mg(esine.cosi(ii +i cos9. sin<)i, ) = f 

X X p X X p X L Z 


xxP 

- I r 
xz 

- I 

xz 

C>iP + 

<^zz 

- I ) 
yy 

(Qj^r + 

Riq) = 


yy<3 

- 21 

xz 

Rj^r + 

(^xx 

“ "zz> 

RlP = 

m 

y 


(11) 

ZZ^ 

-Ip 

xz^ 

^ ^xz 

(Q^r + 


>1 

>1 

H 

+ 

- "xx> 

QlP = 



In equations (10) - (11) , the quantities with subscript "1" refer 
to reference values and the quantities fy» ^z^' ^"'x' ^y' } 
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refer to perturbed values of the externally applied forces and 
moments respectively. 

Equations (10) - (11) represent the uncoupled equations of 
perturbed motion in body axes of an elastic airplane consistent with 
the assiamptions necessary to achieve decoupling. They are to be 
considered further in another section. 

Elastic Body Motion 

In addition to the six equations for the rigid body motion for 
the body axis system, there exist 3n equations of motion which reflect 
the internal equilibrium of the n lumped masses making up the total 
aircraft mass. Here, it is convenient to introduce influence coeffi- 
cients into the analysis which are obtained from a solution ot the Navier 
equations of elasticity. This requires that we consider the deformations 
d| measured relative to a coordinate system (x, y, z) with origin fixed 
to a point coincident with the location of the aircraft center of 
mass before deformation (see figure 2) . After defoinnation, this point 
translates a distance given by the vector ^ relative to the center of 
mass which is assumed to have its position unaltered by elastic defor- 
mation. The coordinate system fixed to point also rotates through 
an angle which can be represented by a rotation vector relative to 
the initial orientation of the coordinate system before deformation 
or loading. Then, the deformation d^ measured relative to a co- 
ordinate system X Y Z with origin fixed at the center of mass of 
the deformed aircraft P, consists of the deformation d^ measured in the 
X, y, z system in addition to the translation ^ and rotation ^ x r^. 
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Figure 2. 



(prior to deformation) 


- Vector Diagram illustrating relationship of to d| 
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That is 


d. = 
— 1 


d! + d + 0^ X r. 
— X — o -o — 1 


(x~*lf2f • • 


( 12 ) 


where is the location of the ith particle relative to the center 
of mass before deformation. ' i 

The displacement field d| is linearly related to the applied force 
field (consisting of inertial, gravitational, and aerodynamic forces) 
for elastic deformations which satisfy a linear stress-strain relation- 
ship (reference 1) . That is, 
n 

= I C. . • Q. , (i=l,2,...n) (13) 

1 j=l 3 


where n is the niimber of lumped masses and C . . 

^ J 

coefficients associated with the mass pair m^, 

C are expressed according to 
11 


is a dyadic of structural 
m^ . The components of 


C . = i IC 

Ij — ‘ xxij 


i + C . . i + C . . k] 
— xyij ^ xzxj — 


+ j [C . . i + C 
-l^^yxi: - yyij 


j + C . . k] 
J- yzxj 


(14) 


+ kIC ..i + C ..j+C ..k] 

— zxxj — zyx 3 ^ zzij — 


The physical meaning of a typical coefficient C . . is that it 

xy 1 J 

represents the component of displacement of the ith lumped mass in 
the X direction due to a unit force in the y direction at the jth 
lumped mass. 
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The result expressed by equation (13) can be rewritten as 


{d’} = [ci{Q}, where 


{d* 

= Id;,, d;,, 


'^x2' ^^y2' 

^z2' • 

(Q) 

~ ^^xl' ^yl' 

Qzi' Q 

x2' °y2' 

°z2' • • • 



^11 

^12 • • 

• ^in' 



^21 

^22 ' • 

• C- 
2n 


[C] = 

^nl 

L 

^n2 ' • 

. C 

nn 

where 


C 

XXX j 

C 

xyi 3 

C 

XZl] 


C . . = 
13 

C 

yxi: 

C 

yyi3 

C 

yzi 3 



C 

ZXX] 

c 

zyi3 

C 

ZZl] 


(15) 


Substitution for d| from equation (13) into equation (12) 
results in the elastic displacement field measured relative to the 
center of mass of the aircraft 


d. 
— 1 


d + e 
— o — o 


n 

X r . + I 

j=l 


C. . 
13 


^3 


f ( 1 — 1 r 2 f • • • n) 


(16) 


where 


Qi = 


-m. ( 


d r 
—1 

dt^ 


- a) + Zi ' (i=if2. 


■ n) 


(17) 
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In equations (16) - (17) , r| is the position vector relative to an 

inertial system (figure 1) of the ith mass particle after deformation 

has occurred and is the aerodynamic force acting on this particle. 

In order to develop a system of differential equations from 

d^ : 

equation (16) we must express the acceleration — i in terms of rigid 

dt^ 


body quantities (center of mass motion) as well as elastic displace- 
ments. Since rl = rl + r. + d., where rl is the position vector to the 

— 1 — u — 1 — 1 — 0 

center of mass relative to the inertial system, we have 


dr 


dt 


— X 

T 


dt^ 


dt^ 


+ li) 


dV 


+ d d(r. + d.) 


dt dt^ dt 


^ + dff^i-^ ^ <£i + ^i> ^ 


which further reduces to 


dr;. 

i = V + OJXV +d. + 2o)xd. 

^t2 -k: - -c -X - -X 


+ £ X (r^ + d^) + ^ X [£ X (r^ + d^^) ] 


(18) 


In equation (18 ) , the accelerations u. x d^ and ^ x x (r^+d^) ] 
are assumed small compared to the other terms which maJte up the absolute 
acceleration of particle m^. Another simplification is realized by 
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computing the acceleration ^ x (r^ + d^) as u x which utilizes 
the small elastic deformation assumption. These simplifications 
allow equation (18) to be written as 

>i2 

dr. . _ 

i- v-^ + toxV+uxr. +d. (19) 

2 —c — — c — —X —X ' 

at 

which results in a value of approximated as 


Q. F. + m.g - m.IV + (oxV+a)xr.+d.] (20) 

—X - —X x^ X — c — — c — — X —X 

The governing equation for the elastic displacement field is then 


d . = d + e X r . 
— X — o -o — X 


+ I C. . • IF . + m.g_ - m. (V + uxV+uxr. +d. )] (21) 

J J J J J 


To be consistent with the previous rigid body results (equations 
(10) - (11) ) we now will write the equations in their perturbed form; 
that is, we let 


Qi 

= Qii 

+ 

-pi 

F . 

— X 

= ^lii 

+ 


d. 

— X 

= lii 

+ 

d . 
-px 

V 
— c 

= ^1 

+ 

V 

-cp 

to = 

(Oj^ + 

(0 

-E 

) 


14 


2 = 2.1 + 2p 

9 = 01 + 6p 

4i = 

il' = >/'l + 4>p (22) 


where the subscripts "1" and "p" denote steady-state or nominal values 
and perturbed values, respectively. Isolation of nominal and perturbed 
values in equation (20) results in the perturbed force 


Q.=F.+m.g -m. - m . to x V 

X ^p X — cp X 


- m^ X V^p - m^ - m^ (Op X r^ , (i=l,2,...n) (23) 


This result is further expanded into component form through the use 
of the relations 


V , = u, i + V, j + W, k 

—cl 1— 1— 1— 

= P^_i + 

£i = ' (i=l,2,...n) (24) 

V = ui + vj + wk 

u = pi + qi + rk 
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= gle cose, i + {(j) cose,cos(j)T - e sin0-sin((iT)i 

P p ± X p X X 

-t0pSin0j^cos<|ij^ + (jipSinij>j^cos0j^) k] (25) 


where i' k. unit base vectors in the principal body axis system. 
Substitution from equation (24) into equation (23) results in 


Q ■ 


F . 


• • 

d . 


1 0 0 0 2 j “V ■ 

pix 


pix 


pxx 


^i JX 

Q • 

_ 

F . 

- m. 

d . 

- ni. 

0 1 0 -z . 0 X. 

piy 


pxy 

X 

pxy 

1 

X X 

Q • 






0 0 1 y . -X. 0 

piz 


pxz 


pxz 


■'x X 


0 


- m 


i 


R 


1 




u 

V 

w 

p 

q 

r 



0 

0 

0 

0 

-gcos0^ 

0 

m . 

X 

0 

0 

0 

-gcose j^cos^j^ 

gsine j^sin^i^ 

0 


0 

0 

0 

gsiniji Jocose 

gsinO j^cosi^ 

0 


O 

O 

o 


p 




( lf2f . . .n) 


(26) 
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This result can be expressed compactly as 


# 

{Qp} = {Fp} - rm.J{dp} - I'rnjm({V^p} + [Mj^HV^p} + lM2l{r^p}) 

( 27 ) 


where 


(Qp) = 

^®plx' 

Qpiy' 

®piz" 

°p2x 

(Fp) = 

I^plx' 

^ply' 

^plz ' 

^p2x 

(dp) = 

t^plx' 


•^plz' 

*^p2x 


= lu, V, 

' w, P, 

- q» rj 

T 


= ro, 0 , 


V" 


' ®pny' ®pnz^ 


' ^pny' ^pnz^ 


' Vz^ 


( 28 ) 

( 29 ) 

( 30 ) 

( 31 ) 

( 32 ) 


['mj = 


m. 


m. 


m. 


o 




m. 




o 


m 


n 


m 


n 


in_ 


n 


( 33 ) 
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0 0 

o 

-gcose 

0 

0 0 

-gcosej^coS()>j^ 

gsine j^sin(|) 

0 

0 0 

gsin<(ij^cos0j^ 

gsinej^cos(^ 

O 

0 0 

0 

0 

0 

0 o 

O 

0 

0 

0 0 

0 

0 

0 


(36) 


The matrices [* m J and [iji] are referred to as the diagonal mass matrix 

and rigid body mode shape matrix, respectively. 

Further, if we combine the vectors d and e into a vector 

— o — o 

= f'^xo' *^yo' *^zo' ®xo' ®yo' express the elastic 

displacement field from equation (16) as 


{d} =[,()] B + IC]{Q} 

where {d} = [d^^^, d^^, d^^, d^^' V' 
and 


{Q> - [Q^i, QyjL' Qzi' Qx 2' “■ °zn^ 


In terms of perturbed quantities, we have 


(37a) 


{dp} = [4'] {Bp} + [C]{Qp} (37b) 

where {B } represents the perturbed value of {B} from the reference 
P 

condition. 
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The vector can be eliminated from the analysis by assuming 

the mean body axis system coincides with the principal axes system. 

According to Milne (references 2 , 3), the mean axis system is 
defined by the equations 


n 

y m. r . X d . = 0 
i=l ^ 


and 

n 

y m. d == 0 
i=l ^ 


(38) 


(39) 


or equivalently 


J {d} = 0 


(40) 


Then, the perturbation in elastic displacement also satisfies this 
relationship or 

I^]'^['mJ{dp} = 0 (41) 

Equations (38) - (39) state that the deformation motion has zero 
linear and angular momentum relative to the body axis system. 

Ashley (reference 4) interprets the assumption of equations (38) - 
(39) as a statement that the overall motion of the body is described 
by a set of principal axes for the deforming body. 

If both sides of equation (37b) are premultiplied by the matrix 

— T 

[ifi] ['mj, the result is 
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J 


I 


[ir]'^rmj{d„} = UHB } + J [C] {Q } = 0 

P r r 

from which {B } can be solved for as 
P 

{Bp} = -[M]"^[ir]’^rin J [C] {Qp} (42) 

where 

IM] = rin'^I'mJI?] (43) 

Substitution of from equation (42) into equation (37b) gives 

{dp} = [C]{Qp} (44) 

where 

[C] = (I - [?] [M]"^[?]'^['mJ) [C] (45) 


If {Qp} from equation (27) is substituted into equation (44) , 
we find that the perturbation elastic displacement field satisfies 
the differential equation 


{dp} = lC]{Fp} - rcjrmj{dp} 

- IC] I'm.] U] ({V^p} + IMiHV^p} + lM 2 ]{r^;^p}) (46) 


It is possible to express the equation for I dp} in terms of the 
stiffness matrix = ICJ ^ by substitution of iQp) from equation 
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(27) into the expression for t<3p} as expressed by equation (37b) . 

This results in 

(d } = [?]{B } + IC]{F } - [C]I-mJ{d } 

r H r' H 

- [cj rmj ({V^p} + + [M2J{r^p}) (47) 

Again, we wish to eliminate {B } from equation (47) . This can 

P 

be accomplished by the following operations. First, premultiply 
both sides of equatioh (47) by = fC] ^ to get 

[Kii]{dp} = [K^3^]U]{Bp} + {Fp} - rmj{dp} 

- rmj 171 ({V^p} + EMj^]{V^p) + lM2]{r;p}) (48) 

Note that the rigid body perturbation equations expressed by equations 
(10) and (11) can be written as 

lM]({Vcp> + f^l^f^cp^ f”2^^''op^^ = I^]'^{Fp} (49) 

where 

[?]'^{F } = [f , f , f , m , m , m 
lYJ x' y z X y z 

Also note that [iji] J {d } = 0 since Ii)> ] ["m^J {d } = 0. 

P P 
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These two results reduce equation (48) to the form 


from which {B } is solved for as 
P 

{Bp} = (50) 

Substitution of this value for (Bp), which is equivalent to the value 
expressed by equation (42), into equation (48) results in the 
equation 

rmOCd } + [K]{d_} = {F„} 
p P P 

- rm.]lT]({V^p} + [M^HV^p} + [M2Hr-p}) 

(51) 


where 


IK] = [K^^j - U] ( U])"^ 


(52) 


Equation (51) represents 3n scalar equations of motion for the 
n lumped masses. Only 3n - 6 of them are independent; however, since 
only 3n - 6 components of {dp} are linearly independent due to the con- 
straint equation (40) defining mean axes. These 3n - 6 independent 
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equations for internal motion together with the rigid body or body 
axes motion from equation (49) make up the necessary 3n equations 
of motion to describe the dynamics of the elastic airplane. 

The results represented in equations (49) and (51) are very 
difficult to treat mathematically without further simplification. 
Before preceding to the various approximate forms for these equations, 
it is of interest to consider them in more detail. To demonstrate 
the roles of the stability derivatives associated with rigid versus 
elastic airplanes, we now develop the right hand side of equation 
(49) . The perturbed aerodynamic force associated with the n lumped 
masses is generally assumed to be a linear function of the velocity 
and acceleration of the aircraft center of mass and of the elastic 
displacement field and its velocity and accelerations fields, i.e., 
from reference 5, 

(Fp) = [A^]{V^p> + [A 2 ]{V^P> + [A 3 ] {dp} + [A^Hdp} + [A^Hdp} (53) 

The matrices [A^] , ... [A^] are referred to as aerodynamic influence 
coefficients which contribute to the stability derivatives. Sub- 
stitution of [Fp} from equation (53) into equation (49) result in the 
following aerodynamic force and moments due to rigid and elastic 
body motions 
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[?]'^lA^]{V^p}, [?]’’lA2]{V^p}, Uj'^[A3]{dp}, 

m'^IA^Hdp}, m'^lAgHdp} 


The nature of the elements of these quantities is demonstrated in the 
following example 


[«l'^[A3^]{V^p} 








1 

1 



• 

1 

0 

0 1 

1 

0 

0 * • • • 
1 

1 

0 

0 

0 

1 

o • 

o 

1 

0 1 • • • 
1 

1 

0 

1 

0 

0 

0 

1 

._Lj_ 

0 

o 

.. i . j .:::. 

1 

-h 

0 

0 

1 

0 

-^1 

~ 1 
^11 

0 

-='2 

^ 21 *” 

1 

1 

0 

-z 

n 

^n 

==1 

0 

-Xi 1 

==2 

o 

- X 2 1 • • • 

1 

z 

n 

0 

^n 

^1 

^1 

° 1 

■^2 

^2 

0 1 * * * 

1 

-^n 

X 

n 

0 










* 


^llxu 
a, , 

^llxv ^llxw ^llxp ^llxq ^llxr 
a__ 

llyu 
a, ^ • 

llyv 

llzu 


a, « 

a, « •••••••••«•••••••••••••• 

12xu 
a, ^ 

12xv 

12yu 


12zu 




a. 


Inxu 
a, • 

Inxr 

Inyu 

a, • 


Inzu 

Inzr 


u 

V 

w 

p 

q 

r 


( 54 ) 
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where the first subscript denotes an element of [A^] , the 
second subscript denotes the element of mass, the third subscript 
denotes the component of force and the fourth subscript denotes the 



or, 



( 55 ) 
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where, for instance, 


f - l/2pSH^(C^ + 2 
u u 

M = l/2pSU^c(C + 2 


(56) 

(57) 


The quantities f and M can be interpreted as the force in the x 
u u 

direction per unit change in velocity u, and moment about the x axis per 
unit change in u, respectively. Also, as a second example, consider 


I*]'^[A3]{dp} = 


^311xx ^311xy ^311xz ^312xx" ' ‘ ^31nxz 


^311yx ^311yy ^31nyz 

^311zx ^31nzz 


[<^] 


T 


^321xx 

^321yz 

^321zx 


32nxz 




^3nlxx 

^3nlyx 

a 

3nlzx 


^3nnxz 

^3nnyz 

^3nnzz 


( 58 ) 
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The subscripts on the elements of lA^] have the following meaning; 

the first subscript indicates an element of [A^J t the second and third 

indicates a mass element pair, the fourth indicates the component of 

force and the fifth denotes the component of displacement. For instance, 

^321zy ^®P^^sents the force acting on panel 2 in the z. direction due to 

a unit displacement of panel 1 in the y direction. We will defer discussion 

— T 

of the stability derivatives contained in the term [(ji] [A_] {d } until 

P 

we have introduced the modal substitution concept. 


Quasi-Static Approximation 

The quasi-static approximation considers only those aeroelastic 
forces dependent upon {<3p} but not upon the time rate of change of this 
displacement field. For example, the inertial forces I^m^]{dp} and 
[^J^[AgJ{dp} and the aerodynamic damping force [^]^[A^j{dp} are con- 
sidered to be equal to zero. This assumption allows for the elastic 
displacement field to be eliminated from the analysis by solving 
equation (46) for (dp) as 

{dp} = lC]{Fp}-[C] rmj It] ({V^p} + [Mi]{V^p> + lM2]{r^p}) (59) 

Substitution of I dp} from equation (59) into equation (53) results in 
{Fp} given as 
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{Fp} = (I - [A3]IC])“^ nA^]{V^p} + [A2]{V^P> 


- [A 3 ] IC] I'm J U] ({V^p} + [Mj^]{V^p> + [M 2 ]{r^p})} (60) 

Substitution for the aerodynamic force given by equation (60) into 

the right hand side of the rigid body equations of motion, equation (49), 

results in the quasi-static form for the rigid body motion as 

[M']{{V^p> + fMi]{V^p} + 1^2^ f^'niA^HV^p} + lA2HV^p}) (61) 

where 

[M'l = [?]'^{[I] + (II] - [A 3 ]IC])"^ [A 3 ] [C] } I'm.] IF] (62) 

U'] = [?]^([I] - IA 3 ]IC]) (63) 

Comparison of equation (61) with equation (49) shows that in the quasi- 
static approximation, the effect of elasticity is to modify the mass 
matrix ["m-.] and the rigid body mode shape matrix, ^ by amounts depending 

upon the elastic properties of the aircraft. 

Modal Substitution 

The method of modal substitution is a means of reducing the large 
niimbers of equations of motion associated with the "exact" formulation of 
equations (49) and (51) while at the same time, uncoupling the degrees 
of freedom of the elastic motion. To this end, we introduce the trans- 
formation 

{dp} = {(|)(x,y,z)} u(t) (64) 
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where {<ti(x,y,z)} is the free vibration mode matrix of constants of 
dimension 3n x 1/and u(t) is the generalized elastic scalar displace- 
ment associated with this mode shape. The values of {(ji} are determined 
from the eigenvalue problem associated with the invacuum motion of the 
aircraft; that is, the deformation of the aircraft vibrating in the 
absence of external forces. The values of u, however, are evaluated 
from the full equation of elastic motion expressed by equation (51) . 

The elastic displacement field for invacuxam motion, neglecting 
structural damping, is governed by the equation 

[■'mj{d } + [K]{d } = 0 (65) 

b' Mr 

In terms of { ij) } and u, equation (65) for each mode is written as 

^rmj{<^} = - [K]{$} 

2 

which requires that ^ be set to a constant say -X , then 

IK] {if.} = X^I'mJ {<})} (66) 

2 

Associated with equation (66) are 3n-6 distinct eigenvalues X^ 
and corresponding eigenvectors The elastic displacement field 

is assumed to be a linear combination of these modes as 

3n-6 

{d } = I {4'.}u. (67) 

P i=i 1 ^ 
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where is the amplitude of the mode shape This is further 

written as 


{dp} = U]{u> 


( 68 ) 


where 


[ 4 ’] — ^ '^2 ^ ^ ^ '^3n~6 ^ ^ 

{u> = '^2'--*"3n-6^'^ 


Hence [cji] has dimension (3n) x (3n-6) and {u} has dimension (3n-6)xl. 

Since the mode shapes correspond to free vibration, it is found 
that the eigenvectors are orthogonal as indicated by the equations 


{4>j}'^rm^]{4.^} = 


(69) 


or 


14}'^ (K] [ 4 ] = ['KJ 
[4]''’['mv] [(()] = I'in.] 


(70) 


The matrices ['•m^] and ['Kj are the generalized mass and stiffness 
matrices, respectively. 
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Substitution for {d }from equation (68) into equation (51) results 
in the expression 


rm,] l4.Hu} + [in HHu} = {F } 


(71) 


- HmH U] ({V^p} + + [M2Hr^p}) 


T 


If both sides of this expression are multiplied by [4>J , it follows that 


[<^]'^hm.H 4 ]{u} + [4]^[K] [4>]{u} = [4.]'^{F } 


-[<^] HmJ[<^]({V^p} + [Mj^HV^p} + [M2 J(r^p}) 


(72) 


op' 


The last term on the right-hand side of equation (72) vanishes due to 
the definition of the mean axis. This reduces the equation of motion 
for the generalized displacements to the form 


Pm,]{u} + ['K.Hu} = [ 4 ] {Fp} 


(73) 


where 


{Fp} = [A^ltV^p} + lA2HV^p} + [A3][4Hu} 


(74) 


+ [A^][ 4 i]{u} + [AgJHJiu) 
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The term [A^] I<(i]{u}, occurring in equation (74) is now considered 
to show the form of some of the additional stability derivatives due to 
elastic effects. The force [i|i] [A^JEifiJCu} discussed in equation (58) 
in terms of the physical displacements is now written in terms of 

generalized displacements as 


[*]'^[A3] M {u} = 


where, for example, 

l/2pU,^S/-C , (i=l,2, . . .3n-6) 

^ u . 

X 
C 

1 

To summarize the results of this section, we outline the sequence 
of computations to arrive at the elastic displacements (dp): 




33 



1) Compute [ <t>] , the mode shape matrix, from the 3n-6 eigen- 
functions of the relation 

[K] = X^['mJ 

If is the ith eigenfunction, then 

[(|i] = [ {<(i^ } { i|>2 } • • • 

2) Compute the generalized mass and stiffness matrices ['•m^] 
and ['"Kj from the orthogonality equations 

U] = ['KJ 

3) Compute the generalized elastic displacement vector {u} from 
the equation 

hm^]{u} + hK^]{u} = [<p]'^{F } 

4) Compute the physical elastic displacement vector {^3^} from 
the equation 

{dp} = [(^]{u} 

Residual Flexibility 

In this formulation, the free vibration mode shapes matrix and 
amplitudes are separated into dynamically retained modes and dy- 

namically deleted modes Iij> 2 J • This means we take the modal substi- 
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I 


tution formulation and omit terms or forces such as ^^^ 2 ^ ' ^^4^* 

[<(12] {^2}? and [Ag] [i))2J {U2 } j where ['•m2^j the generalized mass matrix 
associated with the deleted modes. Here, we have treated the deleted 
modes as we did the total displacement vector in the quasi-static 
formulation. 

We write {dp} as 

{dp} = [(fij^Huj^} + [412HU2} ( 75 ) 


where 


[<!>] = [[■). 

ll/ U'2^^ 

; {u} = 

{Uj^} 

{U 2 } 

^ • 



Further, let 








10] 



[0] 

['K,} = 

[0] 


/ C^in J = 

[0] 

rS2.] 


“ 




_ 


Then, we have (substituting the above relations into equation ( 73 )) 



{u^ } 


to] 


I 

rH 

■e- 




.{u2 3 


[0] t'K2J"^ 


t^23^ 


. - 




. 


{{Fp} - ['m,] [.(.j^Huj^} - ['mxl [4-2] {^2}} 
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Solve equation (76) for {^ 2 ) to get 


-1 


{U2> = -[•'K 2 ,] 


-1 


( 77 ) 


-['K2^1 [<^ 2 ]’^ I'm-] 1^2^ 


This equation is simplified by noting the following 
[<^2]'^['ni.] = 0 

- I'K2^] I'(>2]'^I''m.^j [<^ 2 ] {U2> = 0 

The latter result follows from the basic assumption of residual flexibility. 
This results in the approximation 

{U2> S I'K2^]"^U2^^^^P^ 


where 


(Fp)? lA^HV^p) + lAjHVpp) + 


+ [A^] [<( 12 ] 1^2^ [A^] U^] {u^} + [Ag] {Uj^} 


(79) 


Substitution of {U 2 } from equation (78) into equation (79) results in 
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{Fp} = lA^]{V^p} + [A2J{V^P> + [A3 U^^Hu^} 

+ [A^] {u^} + [Ag] {Uj^} (80) 

+ IA3] I().2] PK2 J ['*2] 

The last term in equation (80) which explicitly relates to the deleted 
mode shape can be rewritten in terms of the identity 


[,fr2H''K2^]"^I<).2]'^{Fp} = [Cj^lKFp} - ["m J {u^} 

- IC] r-m J U] ({V^p} + tM^HV^p} + fM2]{r;p}) 


(81) 


where 


[Cj^] = [C] - [<(. 3 ] (82) 

This identity is obtained by equating the expression for f^p} from 
equation (46) to the expression obtained by using equation (76) . 

Substitution of equation (81) into equation (80) results in 
the following expression for (Fp). 
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{Fp} = ([I] - IA3] [Cj^])"^[tA^HV^p} + [A2J{V^p} + [A3] {U^} 

- [A3] ([C] I'm^] [?] {{V^p} + [M3^]{V^P> + [M2]{r^p}) 

(83) 

+ [Cj^] rm.] [*^]{u^}) + [A^] [,(. 3 ^] {u^} + [A 5 ] [^^] {Uj^}] 

This form for the perturbed aerodynamic force results in the following 
equations of motion for the residual flexibility approximation to the 
elastic airplane: 

Rigid Body Motion; 

[MKtv^pl + + IMaHrip}) - 

- tAjUC^n-^UAiHVpp) + lAjHVpp) 

+ [A 3 ] [<f^] {u^} - [A 3 ] [Cj^J [-m J [<f^] {u^} (84) 

- [A3] [C] rm.] U] ({V^p} + [M3^]{V^P> + [M2Hr^p}) 

+ [A^] [ 4 , 3 ^] {Uj^} + [A 3 J [43^J{U3^}} 
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Elastic Displacement: 


r~mj^ {u^} + {Uj^} = 

[^l]'"([I] - [A3][C^])-l{[Aj{V^p} + [A2]{V^P> 

+ [A 3 ] [<^^nu^} - [A 3 ] [Cj^] [-m J [<j)^]{u^} (85) 

- [A3][C]r'rnjm 

+ [A^] [^^] {U3_} + [A5] [<^j^I {Uj^}} 

Modal Triincation 

The modal truncation form of the equations of motion is obtained 
by representing the modal matrix in terms of only r mode shapes where 
r <3n-6. That is, we let {d } = ['j',]{u^} where 

XX 

[•fij^] = ■[‘f'2^ ••• 

T 

{Uj^} = [u^^, U 2 , ... u^] 

This results in the following equations of motion: 
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Rigid Body Motion: 


[M]({V^p} + [M^HV^p} + [M2l{r;p}) = 


op 


[^]^([A3^]{V } + [A 2 Hv } + [A3] {u^} 


+ [A4] [<^3] {U3} + [A3] [<|.3]{U3}) 


Elastic Displacement: 


[~m^J{u^} + pK^--,] {U3} = 

[<^]^([A 3 ]{Vpp} + [A 2 ]{Vpp} + [A 3 ] [< 1 . 3 ] {u^} 

+ [A4] [(j)^] {U3} + [A3] [<|.3]{G3}) 

PARAMETER IDENTIFICATION 

The parameter extraction process will be discussed for the simpler 
case of an elastic aircraft perturbed from symmetric, steady, level, non- 
accelerating flight. In addition, the equations of motion used to re- 
present the aircraft dynamics will be the modal truncation formulation 
where in vacuum modes are retained. For this case, the rigid body motion 
or mean axis motion from equation (84) is given as 
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Mu + Mge = f u + f w+f 0 + f u 

P P 


c 




m 

+ I 

i=l 


(f 


u . 

1 


+ f 




u. + f 
1 X 


u . 




( 86 ) 


-MU, 0 + Mw 

1 P 


U + f w 
z z 

u w 



+ f 



+ f.w + f.0 + y f 6. 

Z Z D Z 1 

w q ^ i=l <5^ 

m 

+y (f u. +f u.+f„ u.) 

x=l u. u. u. 


I 0 =Mu+MW+M0 + M'U + M’W + M* 0 
XX p u w qp u w qp 


m 


+ T M^ 6. + y (M u. + M* U.+ M" u.) 
. 5.x . u. X u. X u. X 

X=1 X X=1 X XX 


and the elastic motion satisfies 


Tn.U. + K.u. = f. u + f. w + f. 0 + f. u 

^‘X X XX X X. xp X* 

u w q u 


+ f.,w + f. 0 + y f. 6. 

H P j = l => 

m 

+y (f. u. +f.. u. +f.„u.) , 

j=l ^U. 3 ^u. 3 ^u. 


(87) 


. .m) 
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In equations (86) - (87) the forces and moments due to control surfaces 
are included. 

The equations (86) - (87) are in dimensional form. For purposes 
of estimation, they are put in nondimensional form as follows; 


Mu + Mge = l/2pSU^^[(C^ + 2 C^^) u/U^ + S 

^ u a 1 


c5 -• *— 0c 

+ c ^ + c , + c + c — ^ 

V 2U J X* ? X* 2 

q 1 ^u 2U, a 2U, q 4U, 


m 


2u. cu. 

+ I c «. * I (C u . ^ + C ^)1 

1=1 0 . 1=1 u. 'c u. 1 u. 2U, 

1 1 1 1 1 


-MU 6 + Mi . 1/2PSU 2 [ (C^ + 2c ) g + g 

^ u 1 a 1 


C0 

+ c + c 

^q ''l 


C + C 

2 z- 2 z 


52 0 


'u 2U-, a 2U^ 


q 4U^ 


( 88 ) 


° ^ 2u. U. u.c 

ur+'^z- r77>) 

1=1 6. 1=1 U. c u. 1 u. 2U-, 

1 1 111 
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Ixx“p - 1/2 ps5Di^KC„ + 2 + C„ 5- 

U a 1 


0_c 


•' -2 

0 c 


+ C + C + C + C 5- 

”*q ^^1 "*u 2U^ "*0 2U^ ’^q 4U^ 


c m _ • - 

2u. u. u.c 

+ I C 6 . + y (C„ — - + C . + C .. — ] 

j ill nil. 1 jiii m.. - m.. U, m 2 ' 


• •! Ill r. ^ • «i III 

1=1 6. 1=1 u. c 

1 1 


“1 2Xi^ 


m.u. + K.u. = l/2pU, S[C + C 

11 11 '1 u. U, u. U, 

1 , 1 1 1 
u a 


+ C 


0 C 

+ C 


\ 2Ui 


"u 2^1 


+ C ^ + C 
2 u, . 2 u. 


6pC 


(89) 


i- 2U, 
a 1 


i- 4U, 

q 1 


m 


2u . u . 

+ .I_ +.I_ (C„ -^ + C.. + C. 


■n — T ^4 •! — 1 ^4 ^4 m •• OTT 

3 1 ig. 3 1 1^ C 1^ 1 1 2U^ 

D D 1 


u . c 


The stability and control parameters that are to be estimated from flight 
test data are the components of a vector 
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wh0 3T0 (k — I f 2 f w • * c) f 

2 

Th 0 numb 0 r of param 0 t 0 rs list 0 d abov 0 is 3m + me + 3c + 15m + 20. With 
the exception of the control derivatives, unindexed quantities are the 
derivatives associated with rigid airplanes while the indexed quantities 
are the additional derivatives introduced through aeroelastic effects. 

Equations (88) -(89) are now rewritten in matrix notation by 
defining a state vector 


z 


[u,W, 0^ 


0p,U 


l'“2' 






T 


and a control vector 


(90) 


6 = 




(91) 
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The equations of motion then become 


Fz 


Bz^ + 


( 92 ) 


where F, B, and G are matrices of stability and control derivatives 
which are given in Appendix A. In the case of rigid airplanes, all 
state components are available for direct measurement. Additionally, 
linear accelerometers are used to augment the linear speed and angular 
velocities to provide for more efficient parameter extraction. In the 
case of elastic aircraft, additional accelerometer placement is 
desirable since the generalized elastic displacements are not available 
for direct measurement. The analytical form of the acceleration re- 
corded by the ith sensor is 



a 

z 


u 


w 


m 



g 


0 


P 


+ I 

j = l 




U 0 + I 

^ j=l 



(93) 


The form of the sensor measurements, in general, is taken as a linear 
function of the state, that is. 


Z = + _n 


(94) 


where is assumed to be a white noise process of zero mean and co- 
variance R, and H is either the identity matrix or some variation on 
it depending on whether or not the state vector is augmented to include 
the linear accelerations. 
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The maximum likelihood estimate of the parameter vector £ is 
given as (ref. (6)) 

N -N 

£ = £° + [I A^(t^)H^R"^H A(t^)]"^[[ H A(t^) r"^ ^(t ■ ) ] (95) 

i=l i=l 

N 

R = v(t^) v^^(t.) (96) 

In these equations, N is the nvunber of measurements, A is the sensitivity 
9 z 

matrix 2 i(t) the difference y^(t) - ^(p®,t). The asymptotic 

value of the covariance of the error in the estimate is 

N 

E [(£-£)(£-£) '^] = A'^(t^)H'^ R~^ H A(t^)]"^ (97) 

which is the Cramer-Rao lower bound on the covariance of the estimate 
(reference ( 7) ) . 

Several practical problems seem likely in the implementation of a 
computer program to estimate the parameters of elastic aircraft. Even 
for the simple case of symmetric, level flight the number of parameters 
is extensive and is likely to contribute to convergence problems due 
to lack of understanding of the physical significance of the aeroelastic 
derivatives. In addition, considerable computer effort is required to 
provide supporting information in the way of initial estimates on the 
parameters as well as the flexibility matrix and normal mode shapes. 
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CONCLUDING REMARKS 


This effort has presented in a tutorial fashion the manner in 
which the complex phenomena of aeroelastic motion alters the equations 
of motion and stability parameters extraction processes of conventional 
rigid body aircraft. The various degrees of approximation in the 
representation of the motion as provided by the quasi-static, modal 
substitution and residual flexibility have been discussed in some 
detail. The corresponding physical parameters which enter into the 
analysis by each of these approximations are indicated and illustrated 
for the case of symmetric, level, unaccelerated flight. It is anti- 
cipated that potential problem areas exist in the implementation of 
a computer program to extract the large number of stability and control 
derivatives . 
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APPENDIX 


ELEMENTS OF STABILITY AND CONTROL 


The elements of the matrices F, B, and G are 


F: Dimension (4+2m) x (4+2m) 


^11 - ^ ■ 2^^ ^x- 


^12 - - ^ ^x- 


F =0 

13 ^ 

1 —2 

F, . = - ^ c S C 

14 4^ X' 


— 0/ (i— lf2/**.m) 

^1' 4+m+i ^ ‘^x- ' (1=1 f 2,... m) 

^i 


^21 2^® ‘^z* 

u 


Pjj - M - joS o C^, 


F 23 - 0 

F24 - - |cs c^. 


^2^ 4+i ~ ^ > (i— lf2^»«.in) 

^2'4+m+i "" ~ 2^® ^ ^z" ' (i=lF2,...m) 

'^i 

1 —2 
= ^ S C„. 


31 2 


m* 

u 


F32 . - |pS =2 c„. 


MATRICES 
listed below; 
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0 


33 

34 ~ "^xx 


F,. = I.... - S C. 


^3'4+i = 0 


"•q 


(x — Xf2 / • « altl) 


’^3'4+m+i ■ ^m" 

“i 


^i'l == - c 

X. 


(x~'Xf2f • • ain} 

(i=4,5, . . .3+m) 


u 


F.,_ = -^ScC 
1 2 2 u . 

1 • 

a 


^i'3 = 0 

'■i’4 - - i »S Cu. , 


(i=4 , 5 , . . . 3+m) 


(i=4 , 5 , . . . 3+m) 


^'4+.j = 0 


(i=4,5, . . ,3+m; j=l,2,...m) 


^i'4+m+j ®ij ■ 2^° ® 


(i=4 , 5 , - . . 3+m; j=X , 2 , 


Uj 


^4+m,j ~ '^3j' ( j-lr2, . . .4 + 2m) 

^4+m+i,j "" ‘^4+m+i,j ' (i=l,2,...m, j=l , 2 , . . , 4+2m) 

B: Dimension (4+2m) x (4+2m) 


B - i«SU (2 c + C ) 

U 


Bi2 = ^PSU^ 


®13 = ■ Mg 
®14 - = S 


Bw4 


pu^s 


'l'i+4 ^X 

2c u . 

X 


C , (i=l , 2 , . . .m) 


. .m) 
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®l'4+m+i ^x- 


(1=1 , 2 , . . .m) 


'^i 


B21 = |pSU^ ( 2 C^^ -f 


B22 = ""z 


®23 = 0 

B = (M + c )U^ 

, 2 , 


®2'4+i 




2c ^u. 

1 


®2'4+m+i ■ 2^^1® ^z- 


u . 


( X — X f 2 f • • ■ m ) 

(i=l , 2 , . . .m) 


®31 ■ 2»“l S = 


u 


®32 - ° '"m 


B 33 . 0 

B34 - jtSUj c 


B3M+i ■ 


u . 

X 


®3'4+m+i 2 ^^ ^ 


®i'l = 2^^1^ 


®i'2 " *^u. 


(i=l , 2 , . . .in) 

f (i”lf 2 ^... m ) 

(i=4 , 5, . . . 3+m) 
(i=4 ,5, . . .3+m) 


B. ,, = 0 
1 -5 


(i=4 , 5 , . . . 3+m) 


f ^ c ^ f (x— 4 r 5 f • • • 3+ni) 
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pufs 


j+4 


— C - K. 6. (i=4 , 5 , . . . 3+m; j=l,2, 

— U. llj r J f r 


2c 


u . 


r j ^ ^ r 5 r • • • 3+m j 3 — lf2r«««m) 

®4+m,j " '^4j ' (j=l,2, . . .4 + 2m) 


B. , . = 6 . . 
ID ID 


(i=5+m, . . . 4+2iti; j=l , 2 , . . . 4+2m) 


G: Dimension (4+2m) x c 


Gij = , (j=l,2,...c) 

^2 i ~ 2^ ^ (j“l/2f«..c) 

6 . 

D 

G « . ~ U - S c C f ( 3 “lf 2 ,.«.c) 

j D 2 i m . 

j 

'^i+3,j ® ^u. ' (i=l,2,...m; j 

D 


= 1 , 2 , . . . c) 


G. . . = 0 , (i=m,m+l , . . . 2m; j=l,2,...c) 

1 + 4 , J 


. .m) 
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